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We show that the semi-classical model of conduction breaks down if the mean free path of charge 
carriers is smaller than a typical extension of their wavefunction. This situation is realized for 
sufficiently slow charge carriers and leads to a transition from a metallic like to an insulating like 
regime when scattering by defects increases. This explains the unconventional conduction properties 
of quasicrystals and related alloys. The conduction properties of some heavy fermions or polaronic 
systems, where charge carriers are also slow, present a deep analogy. 
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The semi-classical Bloch-Boltzmann theory, which 
plays a fundamental role in our understanding of conduc- 
tion in solids, has limitations such as for example mag- 
netic and electric breakdown l] or quantum interferences 
in the diffusive regime Q that are well known and have 
been intensively studied. The present work focuses on 
another limitation that has received little attention and 
that occurs when charge carriers have sufficiently small 
velocities. 

Indeed the semi-classical theory of conduction in crys- 
tals is based on the concept of a charge carrier wave- 
packet propagating at a velocity V — (l/h)dE n (k)/dk, 
where ll E n (k)" is the dispersion relation for band n and 
wavevector k. The validity of the wave-packet concept 
requires that the extension L wp of the wave-packet of 
the charge carrier is smaller than the distance Vt of 
traveling between two scattering events separated by a 
time r. On the contrary, if Vt < L wp , a condition that 
can be realized by sufficiently slow charge carriers, the 
semi-classical model breaks down. Here we report on a 
quantum theory that allows to treat on the same footing 
the standard regime where the semi-classical approach 
is valid and the regime of slow charge carriers. We find 
that when the scattering time r decreases a transition 
can occur between a metallic regime for Vt > L wp and 
an insulating like regime for Vt < L wp . As an exam- 
ple we consider specifically a complex metallic alloy: the 
a-AlMnSi phase. Ab-initio band structure calculations 
show that the samples of this system, that have been 
studied experimentally, are in the small velocity regime 
Vt < L wp . This explains their unconventional conduc- 
tion properties. The a-AlMnSi phase is structurally re- 
lated to the icosahedral quasicrystalline phase AlMnSi 
and shares many similar conduction properties with other 
icosahedral phases such as AlCuFe and AlPdMn and 
their crystalline approximants 0, 0, @ Thus the present 
work is relevant for these systems too and gives a strong 
insight in the so far unexplained properties of this class 
of materials. The concepts developped here open also a 



new insight in the physics of correlated systems. Indeed 
recent studies of some heavy fermions or polaronic sys- 
tems |E 13 ) where charge carriers are also slow, show 
that their conduction properties present a deep analogy 
with those described here. In particular a transition is 
observed from a metallic like regime at low temperature 
(weak scattering) to an insulating like regime at higher 
temperature (stronger scattering). 

We come now to the treatment of the conductivity of 
independent electrons. According to the Einstein rela- 
tion the conductivity a depends on the diffusivity D(E) 
of electrons of energy E and the density of states n(E) 
(summing the spin up and spin down contribution). We 
assume that n(E) and D(E) vary weakly on the thermal 
energy scale kT, which is justified here. In that case, the 
Einstein formula writes a — e 2 n(E-p)D(E-p), where E-p 
is the chemical potential and e is the electronic charge. 
The temperature dependence of a is due to the variation 
of the diffusivity D(Ep) with temperature. The central 
quantity is thus the diffusivity which is related to quan- 
tum diffusion. 

For a time independent hamiltonian, the quantum dif- 
fusion of states at energy E can be measured through the 
square of their spreading defined as : 



AX 2 (E,t) = ((X(t)-X(0)) 2 ) i 



(1) 



X(t) is the Heisenberg representation of the position of 
the electron along a chosen axis "x" ■ () E means an aver- 
age over all states with energy equal to E. The diffusivity 
is given by 



D(E) 



• lim* 



dAX 2 (E,t) 
' dt • 



(2) 



AX 2 (E,t) can be computed trough the velocity cor- 
relation C{E,t) = (v x (t)V x (0) +V x (0)V x {t)^ E . C{E,t) 

measures the correlation between the velocities at time 
to and to + 1 and is independent of to ■ It can be shown 



2 



600 



400- 



200 




FIG. 1: Square spreading AX 2 of electrons states with 
Fermi energy Ef versus time t, in the cubic approximant a- 
Al69.6Sii3.oMni7.4. AX 2 is the sum of a Boltzmann term, 
AXl = Vpt 2 , and a non Boltzmann term AX^ B . 



that @ 



dAX 2 (E, t) 
dt 



C(E,t')dt', 



(3) 



which allows to compute AX 2 (E,t) in terms of the ve- 
locity correlation function C (E, t) with the initial condi- 
tion AX 2 (E,t = 0) = 0. For a perfect crystal C(E,t) 
is given by C(E,t) = (C(nk,t)} where \nk) are Bloch 
states with energy E n (k) = E and: 

C(nk,t) = 2^|K, m (fc)| 2 cos((£„(fc)-£ m (fc))^). (4) 

rn 

V ntm (k) are matrix elements of the velocity operator V x 
between \nk) and \mk). In the sum the diagonal ele- 
ments (rn = n) are time independent, whereas the other 
terms depend on t. In the perfect crystal AX 2 (E,t), 
which is linearly related to C(E, t) by J3J), is the sum of 
two terms: 



AX 2 (E,t) 



V 2 t 2 



AX 2 B (E,t). 



(5) 



The first term is the ballistic contribution where "V 2 " is 
the average of \V n . n (k)\ 2 over all states \nk) at the en- 
ergy E. The semi-classical theory is equivalent to take 
into account only this first term. The second term in 
iJSJ, AX^ B (E,t), is a non ballistic (non Boltzmann) con- 
tribution. It is due to the non diagonal elements of the 
velocity operator (m ^ n, in and it describes the 
spreading of the wavefunction in the unit cell. It can 
be shown that AX^ B (E,t) < L 2 , where L is the unit 
cell length along the x axis. In presence of a disorder 
characterized by a scattering time r such that the non 
ballistic term dominates i.e. AX^ b (E,t) > V 2 t 2 , the 
semi-classical theory breaks down and the conductor is a 
non standard metal. 

An example of such a crossover from a standard (bal- 
listic) behavior to a non standard (non ballistic) behavior 
is shown on Fig. \T\in the case of the Fermi energy states 
of the complex phase cv-AlMnSi (see below for more de- 
tails) . The non Boltzmann term AX^ B (t) increases very 



rapidly and saturates to a maximum value of the order of 
the square size of the unit cell. At small time t, the stan- 
dard Boltzmann contribution, AX B (t) = V 2 t 2 , is smaller 
due to a very small velocity V of the electrons. Thus a- 
AlMnSi is a non conventional metal at these time scale 
i.e. when the scattering time r is less than ~ 3—4 10 _14 s. 
By contrast in pure Al, AX B is typically 10 3 times larger 
and AX^ B is about 10 times smaller. Therefore the non 
Boltzmann contribution is negligible at these time scales 
and Al is a standard metal. 

For a system with disorder (static disorder and/or tem- 
perature dependent scattering), we compute the diffusiv- 
ity within a relaxation time approximation |9| : 



C'(E,t) = C(E,t) exp(-t/r), 



(6) 



where C is the velocity correlation function of the perfect 
crystal and C is that of the system with defects. A relax- 
ation time r is defined, which meaning is the following. 
After |JBJ| and J3J|, at a time scale t < r the propagation 
of wavefunctions is not affected by disorder and is that 
of the perfect crystal, but at a time scale t > r the prop- 
agation becomes diffusive due to the disorder. ©, © 
and J2J lead to the expressions of the diffusivity D and 
of the dc-conductivity a of the system with defects: 



2. , 



D = V 2 t + 



a = e 2 n(E F )V 2 T + e 2 n(E F 



L\t) 



(7) 
(8) 



Thus D and a have two parts. The first ones, named in 
the following D B and <7b , are the usual Bloch-Boltzmann 
results, and the second ones, named -Dnb and <tnb, are 
due to the non Boltzmann contributions. L 2 (t) is a 
proper average of AX^ B (Ep,t) on a time scale r. The 
second term in the right hand side of © is equivalent to 
the conductivity of Anderson insulators in the Thouless 
regime [ToL 1 1 1| . Indeed in this regime the wavefunction 
spreads between two inelastic scattering events but the 
spreading saturates to the localization length. Here the 
localization length is replaced by the large time limit of 
L(t) which, as can be shown, is bounded by the unit cell 
length L. 

We come now to a detailed analysis of the unusual 
transport properties of the cubic 1/1 icosahedral ap- 
proximant AlMnSi . 1 2| . These results will be compared 
with those of pure Al (f.c.c), cubic Al^Mn ^| and or- 
thorhombic AlgMn 13] crystals, which have a metallic 
behavior. For all those systems, the electronic structure 
calculation have been performed 0, 0, by using the 
self-consistent Tight-Binding Linear Muffin Tin orbital 
(TB-LMTO) method Starting from the computed 

eigenstates \nk), we use equations ©and Q to com- 
pute the velocity correlation function without disorder, 
C(E, t) and the square of the spreading AX 2 (E, t). The 
average ( ) B is obtained by taking the eigenstates of each 

k vector with an energy E n (k) such as: E — AE/2 < 
E n (k) < E + AE/2. AE is an energy resolution in 
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FIG. 2: (A) LMTO total DOS n(E F ) and (B) diffusivity 
D(Ef) in the cubic approximant a-Al69.6Sii3.oMni7.4 (r = 
3 10 -14 s). D is the sum of a Boltzmann term, Db, and a non 
Boltzmann term, Dnb (see equation 0). 



the calculation. The number of k-points per first Bril- 
louin zone included in this average calculation is Nk- For 
a given AD, Nf. is increased until C does not depend 
significantly on N^- Therefore, there is not adjustable 
parameter in our calculation. For Al (f.c.c), Alj^Mn, 
Al 6 Mn and a-Al 69 . 6 Sii 3 . Mni 7 .4, AD is equal to 0.272; 
0.272; 0.272; 0.0272 eV, respectively, and N k is equal to 
1,728,000; 64,000, 18,000; 32,768; respectively. 

For the a-AlMnSi phase, we use the experimental 
atomic structure [T^ and the Si positions proposed by 
E. S. Zijlstra and S. K. Bose 15] to calculate trans- 
port properties in an a-phase with composition: a- 
Al69.6Sii3.oMni7.4. This phase contains 138 atoms in a 
cubic unit cell: 96 Al atoms, 18 Si atoms, and 24 Mn 
atoms. In Fig. |3 the total DOS n of the a-AlMnSi phase 
is presented versus the energy. Its total density of states 
is characterized 14] by a depletion near the Fermi energy 
Dp, called pseudo-gap, which is observed experimentally. 
This is due to a Hume-Rothery mechanism ^| 0] of 
band energy minimization that is strong in quasicrys- 
tals and their approximants [l7]. Following the Hume- 
Rothery condition, it is expected that the most realistic 
value of Ep in the actual a-phase corresponds to the min- 
imum of the pseudo-gap, i.e. Ep — -Ef(lmto) = — 0.163eV 
for our calculation. Within a relaxation time approxima- 
tion the diffusivity D(E, r) is calculated. The Db val- 
ues (Fig. |5Jl are similar in magnitude to those obtained 
by Fujiwara et al. |14| for the idealized approximant a- 
Aln4Mn24 approximant. Dnb is almost independent on 
E, whereas the Db values depend strongly on E and is 
particularly small in the pseudo-gap. 

We present now the results of the conductivity cal- 
culations of the a-AlMnSi phase assuming the value of 
the Fermi energy given above i.e. Ep — -Ef(lmto) — 
— 0.163 eV. The Fig. shows the predicted static con- 
ductivity (dc conductivity) versus the inverse scattering 
time. In the case of a-Al69.6Sii3.oMni7.4 two regimes 



FIG. 3: Ab-initio dc conductivity, a, in cubic approximant 
a-Al69.6Sii3.oMni7.4. a is the sum of a Boltzmann term, <tb, 
and a non Boltzmann term, <tnb (see equation @). Inset: 
Real part, Kea(uj), of the optical conductivity for three r 
values, uj is the pulsation. Simple line, r = r*/3; line with 
triangles, r = r* = 3.07 10 -14 s; line with crosses, r = 3r*. 



appear clearly: the metallic regime (Boltzmann regime) 
at large scattering time, r > r*, and the insulating like 
regime (non Boltzmann regime) at small scattering time, 
r < t* . t* = 3.0710~ 14 s is defined as the time for 
which the Boltzmann and non-Boltzmann contributions 
are equal. As expected from our model, <7nb IS almost 
proportional to 1/t. Therefore, in the non Boltzmann 
regime, the conductivity increases with disorder as ob- 
served experimentally. For realistic r values, r < t* 
[l8| . ctnb dominates and a increases when 1/t increases 
i.e. when defects or temperature increases, a varies 
from 250 (Ocm)- 1 for r = 3.310~ 14 s, to 2000 (Ucm)' 1 
for t = 10~ 15 s. This is consistent with experimen- 
tal results in a-AlMnSi: cr(4 K) ~ 200 (ft cm)" 1 and 
cr(900 K) ~ 2000 (Slcm) -1 and with standard estimates 
for the scattering time in these systems 5] . Furthermore 
for r equals a few 10 -14 s, i.e. when the Boltzmann term 
is negligible, the mean free path is given by the square 
root of the saturation value of AX^ B and is of the or- 
der of 15 A. This is in agreement with estimates in the 
literature As discussed in this means also that 
the systems is far from the Anderson transition despite 
its low conductivity. Within the relaxation time approx- 
imation used here, the optical conductivity is also the 
sum of two terms. One is the Boltzmann contribution 
(diagonal elements of the velocity operator) which gives 
rise to the so-called Drude peak and the other is the non 
Boltzmann conductivity (off diagonal elements of the ve- 
locity operator). As shown in the inset of Fig [21 a Drude 
peak can be identified in the Boltzmann regime, r > r*, 
whereas in the non Boltzmann regime, t < t* , the Drude 
peak disappears. This absence of the Drude peak is ex- 
perimentally demonstrated for several icosahedral phases 
[l!j . We note also that the physics of the non Boltzmann 
regime is similar to that obtained with anomalous diffu- 



4 




le+13 le+14 le+15 

1/T (fT l ) 



FIG. 4: Ab-initio electrical conductivity <j(Ef) versus in- 
verse scattering time 1/r (logarithmic scales): in Al, Ali2Mn, 
ALjMn, and in cubic approximant a-Alg9.6Sii3.oMni7.4. In Al, 
the Boltzmann term dominates, and the model is compatible 
with a simple Drude model (dashed line) le- 
sion in the subdiffusive case 

Finally Fig. 0] exhibits the fundamental difference be- 
tween the cv-AlMnSi phase and a standard metal like Al 
(f.c.c). Ali2Mn and AlgMn phases are somewhat inter- 
mediate and show the first signs of an insulating like 
regime at the strongest values of the scattering rate. 

To conclude, a regime of charge carriers with small ve- 
locities has been identified, where quantum effects play a 
fundamental role and lead to a transition from a metal- 
lic like behavior to an insulating like behavior when the 
scattering rate increases. Quasicrystals and related com- 
plex metallic alloys give an example of this insulating like 



regime. Other systems are known where the charge car- 
riers velocity is small, either in all directions or in specific 
directions, thus we expect that the insulating like regime 
could be observed in other systems. For interacting elec- 
trons, a similar regime should also exist provided that the 
mean free path is smaller than a characteristic extension 
of the quasiparticle wavefunction, which means that the 
Fermi liquid picture breaks down. Several interacting 
electrons systems are good candidates for this small ve- 
locity regime. Recent results suggest this picture for a 
system like CeBg which is a heavy fermion system with 
large mass and thus small velocity. At low temperature 
T the Fermi liquid description is valid but at higher tem- 
perature (T > T*) an incoherent insulating like regime is 
reached in which the resistivity decreases with increasing 
temperature. Within a Dynamical Mean Field Theory 
(DMFT) calculation [24| a disappearance of the Drude 
peak in the incoherent regime i.e. above T* is predicted 
13 . Recent theoretical results for the small polaron prob- 
lem, with narrow polaronic band and thus small velocity, 
have been obtained also with the DMFT. They exhibit 
an analogous behavior Q ■ The static resistivity (dc resis- 
tivity) increases with increasing temperature up to some 
To and then decreases. Again in the incoherent regime 
above T the Drude peak disappears @- A theoretical 
analysis similar to the discussion given here can be made 
also in the DMFT. It should give an additional insight in 
the quantum transport of correlated charge carriers with 
small velocities. 

The computations have been performed at the Service 
Informatiquc Recherche (S.I.R.), Universite de Cergy- 
Pontoise [25|. GTL thanks Y. Costes and J.C. Baccon, 
S.I.R., for computing assistance. 
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